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Abstract: The dynamic responses of finite element model do not 
match with those of the experimental results due to modeling 
and measurement errors. It is impractical to collect the data 
for the full set of DOFs of a dynamic system. And there are 
some cases to have difficulty in measuring the response such 
as slope. This work presents an analytical method to estimate 
the frequency response function (FRF) at the unmeasured 
nodes of the actual system. The method is derived by 
minimizing the discrepancy between the analytical and actual 
FRF data sets dividing into the real and imaginary parts. The 
validity of the proposed method is investigated in a numerical 
application to estimate the unmeasured slope FRF data of a 
beam. 
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1. Introduction 

An accurate dynamic finite element model of a structure is very 
important for structural design and analysis. The frequency 
response of the finite element model does not correspond to 
experimental measurements due to the existence of uncorrected 
noise content and the modeling errors. If a discrepancy between 
the two is found, the analytical results should be modified based 
on the experimental measurements. It is improper to collect a 
full set of dynamic responses, and all response data cannot be 
obtained from vibration testing. The problem can be overcome 
by data reduction or data expansion, i.e., either reduce the 
system matrices to the number of measured DOFs or conversely 
expand the measured response vectors to the full size of the 
finite element model matrices. 

Most of the utilized expansion techniques involve the use of the 
finite element model as a mechanism to complete the 
unmeasured DOFs from the experimental model. Imregun and 
Ewins [1] presented mode shape expansion methods to estimate 
the modal displacements at unmeasured coordinates using the 
incomplete set of measured data together with a corresponding 
finite element model. Le vine -West et al. [2], [3] evaluated the 
robustness and reliability of the Guy an static expansion method, 
the Kidder dynamic method, the Procrustes method, and the 
penalty method along with least-squares minimization 
techniques. Dubbaka and Houghton [4] provided a correction 
method to modify the measured mode shapes by restricting the 
orthonormality and by applying the constraint rules. Based on a 
dynamic modal expansion that minimizes the residual error in 



the eigenvalue equation for each measured mode, Kenneth and 
Francois [5] presented an algorithm for expanding the measured 
mode shapes obtained from modal testing of the full set of 
DOFs of a corresponding finite element model. 
Measured FRF data provide certain advantages over modal data. 
The FRF data contain the effect of the vibration modes and 
dynamic characteristics. A full set of FRF data should be 
expanded using measured data. 

An FRF matrix of the full set of DOFs measured by experiments 
is used to predict parameter matrices of the dynamic system. It 
has been reported [6] that FRF data provide more information 
than modal data, as the latter are extracted from a very limited 
frequency range related to resonance. Friswell and Penny [7] 
proposed an approach to reduce the model order so that the 
stated estimation process reduces to a least squares problem 
based on the FRFs. Utilizing FRF data measured at specific 
positions, with dofs less than that of the system, as constraints to 
describe a damaged system, Rahmatalla et al. [8] identified 
parameter matrices such as mass, stiffness and damping 
matrices of the system, and provided a damage identification 
method from their variations. 

Considering the inconsistency between the analytical and 
experimental responses due to the modeling and measurement 
errors, the FRF data should be corrected for the subsequent 
design and analysis. This work presents an analytical method to 
estimate the FRF at the unmeasured nodes of the actual system. 
The method is derived by minimizing the discrepancy between 
the analytical and actual FRF data sets dividing into the real and 
imaginary parts. The validity of the proposed method is 
investigated in the FRF data expansion to estimate the 
unmeasured slope FRF data of a beam. 

2. Formulation 

Dynamic responses are expressed in the time domain and 
frequency domain. For linear systems there is little loss of 
information going from the time domain to the frequency 
domain. The FRFs in the frequency domain are measured 
instead of displacement and force individually. The dynamic 
responses can be initially obtained from a mathematical model 
such as a finite element method. The dynamic behaviour is 
approximately discretized for n DOFs and can be described by 
the equations of motion as 

Mii + Cu + Ku = f(t) (l) 
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where M and K denote the nxn mass and stiffness 



matrices, 



» = [, 



7 



Ce/T 



is the 



damping matrix, and f(f) is the nx\ excitation vector. 

Inserting U = \]e lQx and f = Fe* nt into Eq. (1) and 
expressing it in the frequency domain, it follows that 
D ( , ( .(Q)U(Q) = F(Q) (2) 
where the dynamic stiffness matrix of the initial dynamic model, 
D ini , is written by 



D. m .(Q) = K-Q 2 M + /QC 



(3) 



where the subscript ini represents the initial finite element 
Q denotes the excitation frequency, 



model, i = V — 1 

F(Q)=[F 1 f 2 



and 



u(n)=[t/! u 2 ■■■ uj represent the Fourier 

transform of the force and response vectors f and U , 
respectively. Equation (2) is valid for an excitation frequency 
Q . The nxn FRF matrix is defined as H in . , 

H^K-Q^ + ZQC]" 1 (4) 

where H ini i - in the FRF matrix denotes the displacement 

response measured at location i due to the unit impulse force 
input at location j. 

The FRF matrix of the initial model is expressed by dividing 
into the real and imaginary parts as 



H w (n) = H WifB (n)+ffl[ Wifcl (n) 



(5) 



The FRF matrix corresponding to a full set of DOFs of the 
dynamic system can be numerically obtained by dividing into 
the real and imaginary parts. 

If variable parameters of mass, stiffness and damping are 
randomly added due to the modeling and measurement errors, 
the dynamic equation of Eq. (2) is modified as 
D mod (Q)u(Q) = F(Q) (6) 

where D mod (Q) denotes the dynamic stiffness matrix of the 

modified system by 
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D mod (Q) = (K + AK)+ /n(C + AC) - Q 2 (M + AM) (7) 

And the FRF matrix of the modified system, H mod , is written 
by 

U(Q) = H mod F(fi) (8) 

where 

H mod (Q) = [(K + AK) + /Q(C + AC) - Q 2 (M + AM)]" 1 

The modified FRF matrix is expressed by the sum of the real 
part and imaginary part as 



H mod (Q) = H mod/e (Q)+/H mod . (n) 



(9) 



The FRF data at unmeasured positions should be estimated by 
expanding the measurement data. Assume the FRF data of the 
displacement responses due to the unit impulse at m{m < n) 
different locations less than the number of total DOFs of the 
system are collected. The mxm measured FRF data set, 
H mod , can be expressed by 

AU = BF (10) 
where A is an mxn Boolean matrix to define the 
displacement measurement locations and B represents the 
mxn FRF matrix including m receptances due to the m 
impulse excitations, U denotes an nxm displacement 

T 

response matrix, B = U , and F is the nxm excitation 
matrix to define the unit impulse at m locations. The matrices B 
and U are also expressed in the real and imaginary terms of FRF 
matrix as 

B = B +/B. , U = U +/U. (11) 

re im ' re im \ 1 1 s 

This work presents an analytical method to expand U or B 
including the experimentally measured FRF data to a full set of 
DOFs. The measured FRF data represent as the constraint 
condition for the expansion. The expansion is performed by 
dividing into the real and imaginary parts, and a least square 
method utilizing cost functions of 

^ |"^7ra,re \^re ^^ini,re )|| |^"zm,z'ra \^z'ra "^"zra,z'ra)|| 

is applied. In order to utilize Eq. (10) into Eq. (12), Eq. (10) is 
modified as 



AH H U 

ini, re ini, re re 



B F 



(13) 
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ATT IT TT Ti TT 

inijm inijm im im 



(14) 



Letting R x = AH^ re and R 2 = AH)' n f M , and solving the 
result with respect to H2j/*U_ and H^'jjXL, , 



respectively, we obtain that 
H^U w =R 1 + B fB F + (l-R 1 + R 1 V 1 

H^U /m =R^ m F+(l-R^R 2 )y 2 



(15) 



(16) 



where the superscript '+' indicates the Moore -Penrose inverse, 
and y { and y 2 are arbitrary matrices. Equations (15) and (16) 
have an infinite number of solutions due to the presence of the 
arbitrary matrices y { and y 2 . The expansion is carried out by 
providing the condition to minimize the cost function of Eq. (12) 
of all solutions of Eqs. (15) and (16). Utilizing the condition to 
minimize Eq. (12) into Eqs. (15) and (16), and solving it with 
respect to the arbitrary matrices, respectively, they are derived 
as 

y, =(l-R^R 1 )(l-R:B re F) + R^R 1 s 1 (17) 



y. (i-r r )(i-r;b,„k)+r r 2 s : 



(18) 



where S l and S 2 are also arbitrary matrices. Utilizing Eqs. 
(17) and (18) into Eqs. (15) and (16), respectively, and arranging 
them, the following equations are obtained as 



H^U w =R 1 4 B w F+^-R 1 + R 1 ) 



(19) 



H:„SU,. m = R^B m F + (i - R;R 2 ) (20) 

Arranging Eqs. (19) and (20), and pre -multiplying both sides of 



the resulting equations by H)^ re and H^- Z - m , respectively, 



i 1/2 . . 

mi,im 



the n x n resulting FRF matrices to be expanded to a full set 
of DOFs can be written by 

U„ = + H Lt ( AH Lt ) + ( B " AH ^ ) ( 21 ) 
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U /m = H)^. m + H)if M (AH^. 2 . m ) (b - AH^. 2 /m ) (22) 

The numerical results of Eqs. (21) and (22) are expressed by the 
real and imaginary parts as 



U re U ^ re re + iU re j m 



(23) 



(24) 



The resulting FRF matrix can be expressed by summing Eqs. 
(23) and (24) as 



^ mo d \^re,re ^ im,re)~^~ r e,im ^ im,im) 



(25) 



Equation (25) represents the n x n FRF matrix expanded from 
the mxm measured FRF data set. It is shown that the 
proposed method can be explicitly utilized in expanding the 
FRF data. 

3. A numerical experiment 

A fixed-end beam, shown in Fig. 1, was considered to 
investigate the validity of the proposed method to expand the 
FRF data. The nodal points and the members for the finite 
element beam model are numbered in the figure. Each node has 
two DOFs of the vertical displacement and slope, and the beam 
is modeled by 49 nodes and 98 DOFs besides the boundary 
conditions. A 5-m long beam was modeled using 50 beam 
elements of 100mm. The elastic modulus of the beam was 
2.0 x 10 MPa , and its unit mass and gross cross-section were 
7,860 x 10 " 6 kg/m and 200mmxlOQmm, respectively. 
The stiffness and mass matrices can be established using the 
mechanical properties of the beam. The damping matrix was 
assumed as the Rayleigh damping to be expressed by the 
stiffness matrix and the proportional constant 0.0001. 
The analytical FRF matrix by the given parameter matrices 
doesn't coincide with the one of the actual beam due to the 
modeling and measurement errors. For this numerical 
experiment, the simulated parameter matrices of the actual beam 
were established as follows: 
M 



mod 



= M W (I 



K mod -K ini (l + a k <j k ) 



'mod 



= c w (i 



+ a cj 



) 



a k and a c 



where I denotes the identity matrix, OC m 
represent the relative magnitudes of the error for the mass, 
stiffness and damping matrices, respectively, and C m , <J k 

and <J C are n x n random number matrices in the range 

[—1, l]. Figure 2 exhibits their variation deviated from the 
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analytical parameter matrices. The error magnitudes included in 
the parameter matrices were established as 



a 



0.01 . 



It's not easy to measure the slope responses of the beam due to 
the external excitation. Using the proposed method, this study 
was performed to estimate the receptance FRF data 
corresponding to all slope DOFs from the measured 
displacement FRF data. Figure 3 represents the receptance FRF 
magnitude due to the impulse excitation at node 5. In the figure, 
the dashed and solid lines indicate the estimated and actual 
value, respectively. As shown in the plots, the displacement and 
estimated slope FRF receptance curves exhibit the similar 
values at the resonance frequency region as the actual values. 
However, the values at the antiresonance frequency represent 
large difference between two data sets. It should be overcome by 
increasing the number of measured DOFs. 




50 @ 100mm 
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Fig. 1 Finite element model of a fixed-end beam 
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Fig. 2 Variation in parameter matrices: (a) stiffness 
corresponding to the displacements, (b) stiffness corresponding 
to the slopes, (c) mass corresponding to the displacements, (d) 
mass corresponding to the slopes, (e) damping corresponding 
to the displacements, (f) damping corresponding to the slopes 
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Fig. 3 Receptance FRF data due to the impulse force at node 5: 
(a) displacement FRF at node 1, (b) slope FRF at node 3, (c) 
slope FRF at node 12, (d) slope FRF at node 11 

4. Conclusion 

Analytically calculated FRF data do not coincide with the actual 
test data due to the construction and measurement errors. And it 
is not practical to collect the FRF data at all DOFs. This work 
presented a data expansion method to estimate the FRF data at 
unmeasured nodes. The displacement and estimated slope FRF 
receptance curves exhibit the similar values at the resonance 
frequency region as the actual values. However, the values at the 
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antiresonance frequency represent large difference between two 
data sets. It should be overcome by increasing the number of 
measured DOFs. 
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